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It is shown how different globally coupled map systems can be analyzed under a common frame- 
work by focusing on the dynamics of their respective global coupling functions. We investigate how 
the functional form of the coupling determines the formation of clusters in a globally coupled map 
system and the resulting periodicity of the global interaction. The allowed distributions of elements 
among periodic clusters is also found to depend on the functional form of the coupling. Through 
f^ , the analogy between globally coupled maps and a single driven map, the clustering behavior of the 

c7^ ' former systems can be characterized. By using this analogy, the dynamics of periodic clusters in 

Cn ' systems displaying a constant global coupling are predicted; and for a particular family of coupling 

functions, it is shown that the stability condition of these clustered states can straightforwardly be 
derived. 
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Q ! I. INTRODUCTION 

u 

rH ■ There has been much interest in the study of the collective dynamics of chaotic systems subjected to global in- 

teractions. Such systems arise naturally in the description of arrays of Josephson junctions, charge density waves, 
multimode lasers, neural dynamics, evolutionary, chemical and social networks |0-p[. The globally coupled map 
(GCM) lattice M constitutes a prototype model for such global-coupling dynamics. It has recently been argued that 
GCM systems yield universal classes of collective phenomena |Q . Specifically, a GCM system can exhibit a variety 
►^ ' of collective behaviors such as clustering (i.e., the formation of differentiated subsets of synchronized elements in the 
\^ [ network) |r|[: non-statistical properties in the fluctuations of the mean field of the ensemble g ; global quasiperiodic 
^^ ■ motion [p|JlC|] ; and different collective phases depending on the parameters of the system ||9| . It has been shown that 
^D a GCM system is closely related to a single map subjected to an external drive and that this analogy may be used to 

t . describe the emergence of clusters in GCM systems in geometrical terms [Q . 

In particular, the phenomenon of clustering is relevant as it can provide a simple mechanism for segregation, 
ordering and onset of differentiation of elements in many physical and biological systems. In addition to GCM 
systems, dynamical clustering has also been found in a globally coupled Rossler oscillators ||l3|, neural networks pSJ, 
C ' and coupled biochemical reactions [H. The interest in this phenomenon has recently grown, since dynamical clusters 
have been observed experimentally in an array of electrochemical oscillators interacting through a global coupling 

-- @- 

^ . In this paper, we investigate the process of cluster formation in general globally coupled map systems by focusing 

on the dynamics of their global coupling functions. In most studies on GCM systems, the mean field of the network 

j^ , has been used as the global coupling function. Here, we study GCM systems subjected to different global coupling 

C^ functions and show how they can be analyzed under a common framework. We investigate how the distribution of 
elements among a few clusters and their periodicities depend on the functional form of the global coupling. 

Section II contains a description of the dynamics of different global coupling functions in GCM systems and a 
calculation of the possible periodicities and cluster sizes when two clusters emerge in these systems. The driven 
map analogy is employed in Sec. Ill to interpret the clustering behavior of GCM systems. In Sec. IV the dynamical 
properties of periodic clusters in systems exhibiting a constant global coupling are predicted; and for a particular family 
of global coupling functions, the stability condition for these clustered states is derived in an Appendix. Conclusions 
are presented in Sec. V. 

II. DYNAMICS OF GLOBAL COUPLING FUNCTIONS 

Consider a general globally coupled map system 

xt+i{i) = (1 - e)/ {xt{t)) + eH{xt{l),xt{2), ..., xt{N)), (1) 
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where xt{i) gives the state of the element i (i = 1,2,..., TV) at discrete time t\ N is the size of the system; e 
is the couphng parameter; f{x) describes the (nonUnear) local dynamics, which in the present article is chosen 
to be the quadratic map f{x) = 1 — rx^\ and Ht {xt{i), Xt{2), . . . , xt{N)) is the global coupling function. We 
shall consider a general class of global coupling functions of N variables such that H{. . . ,xt{i), . . . ,xt{j), ■ ■ ■) — 
H{. . . ,xt{j), ■ ■ ■ ,Xt{i), ■ ■ .), yi,j', that is, H is assumed to be invariant to argument permutations. This property of 
the coupling function ensures that, at any time, each element of the globally coupled system is subjected to the same 
influence of the coupling term. Some examples of coupling functions belonging to this class are 
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H= X =Y[\xt{i)\^. (5) 



The first two examples correspond to forward and backward mean field coupling, respectively, and they have been 
widely used in GCM studies. The third global coupling function is the usual dispersion or mean square deviation of 
N variables, and it may describe systems whose elements do not interact when they are synchronized. This kind of 
global interaction might be relevant in some biological or social systems where the members of a community are driven 
by their deviations from the mean behavior. The last example is the geometric mean. This type of multiplicative 
coupling occurs, for instance, in a system of N sequential amplifiers where the gain of element i is a function of the 
magnitude of its state Xt{i), and H is proportional to the total gain of the system. Many statistical functions of N 
variables share the property of invariance under argument permutations and they could as well be taken as global 
coupling functions in GCM systems given by Eq.([|). 

For some range of its parameters the GCM system in Eq.(n]) reaches an asymptotic collective behavior characterized 
by the segregation of the elements into K clusters, each exhibiting a period P, where the fcth cluster has a number 
Nk of elements, with J2k=i ^k = N. The fraction pk of elements in the kth cluster is pk = Nk/N. The evolution of 
the kth cluster may be described by a variable Xt(fc) which gives the common state of the Nk elements belonging to 
this cluster at time t. The periodic orbit adopted by the state Xt(fc) of the A:th cluster can be expressed as a sequence 
of P values [xi(fc),X2(fc), • • ■ ,Xp(^)]- The specific partition {pi,p2, ■ ■ ■ ,Pk} into K clusters and the specific values 
taken by the periodic orbit of each cluster depend on initial conditions and parameters of the system. 

When a GCM system falls into K periodic clusters, the coupling function H also shows a periodic motion. As an 
illustration of this behavior. Fig. 1(a) shows a typical situation in which a GCM system, with global coupling function 
H = X, displays two clusters, each in period two. In this case, the coupling function follows a period-two motion. 

Collective states consisting of two clusters have recently been observed in an experimental array of globally inter- 
acting chemical elements [|5|. This clustered collective state is interesting to analyze in globally coupled maps. In 
this simple situation, there is a fraction p of elements in one cluster and a fraction (1 — p) in the other cluster. Thus 
the global coupling functions from Eqs.(||)-(||) simplify to: 

H=(x)=pxtil) + il-p)xt{2), (6) 

H = (/)= p /te(2)) + (l-p)/(xt(2) ), (7) 

H = Ax^ 7^(rT^)(^^;^^P^;(I))I, (8) 

H= s=|xt(l)rixt(2)p-^'. (9) 

When the GCM system reaches a two-cluster state, the dynamics of the system reduces to the two coupled maps 

Xt+i(l) = (l-e)/(xt(l))+ei?, (10) 

Xt+i{2) = {l-e)f{xt{2))+eH; (11) 

where H — iJ(xt(l), Xi(2),p) is the reduced, two-cluster coupling function, as in Eqs.(ra)-(p[). If both cluster states 
Xt(l) and Xt(2) fall in period-two orbits, the coupling function H follows, in general, a period-two motion, as shown 



in Fig. 1(a), although H may become constant in some circumstances (see Fig. 1(b) and Sec. IV). Let Hi and H2 
be the values adopted alternatively in time by H in its period- two orbit for a given partition {p, 1 — p}, as indicated 
on Fig. 1(a). The values Hi and H2 depend on i) the functional form of H; ii) the parameters of the GCM system 
Eq. (|l]) ; and iii) the fraction p. 

Consider then several GCM systems with the same parameters r — 1.7 and e — 0.2 but with different coupling 
fmictions such as those in Eqs. (^-(|^). Two clusters in period two can emerge in each of these systems for some 
range of the fraction p. The resulting asymptotic orbits [Hi,H2] of the respective coupling functions are shown in 
Fig. 2(a) as p varies, giving rise to a curve in each case. Note that each function H possesses period-two orbits only 
for a limited range of the fraction p. Figure 2(b) is a magnification of Fig. 2(a) which shows that the dynamics of the 
backward and the mean field coupling functions become equal for p = 0.5; i.e., when the two clusters have equal sizes. 
In this case, both coupling functions reach the constant value Hi — H2 = 0.365. Notice that the dispersion coupling 
function, H = Ax, only displays states with Hi = H2] that is, even when the two clusters in period two may have 
different sizes, this particular global coupling always reaches a stationary value. The curves for the other coupling 
functions are symmetrical with respect to the diagonal in Fig. 2(b), which they cross for p = 0.5. On the diagonal, 
the coupling functions are constant and the two clusters evolve out of phase with respect to each other (Sec. IV) . 

Note also that the different global coupling functions perform a period-two motion only on a restricted region of 
the plane (Hi,H2). It will be shown in Sec. Ill that period-two orbits of any permutable coupling function will fall 
within the dashed contour in Fig. 2. 

In general, a coupling function Hofa GCM system in a collective state of two clusters can reach various asymptotic 
periodic orbits for appropriate initial conditions. Each Fig. 3(a) to 3(d) shows the regions on the space of parameters 
(p, e) for which a coupling function H oi a GCM in a two-cluster state displays different periodic motions. The local 
parameter is fixed at r = 2. Figures 3(a) and 3(b) correspond to the backward and forward mean field coupling, 
respectively. Note the very different distributions of periodic regions for the coupling functions in Figs. 3(a) and 3(b). 
It should be noticed that, besides the collective periodic states for two clusters shown in Figs. 3(a)-3(d), there can 
exist other states in a GCM system consisting of more than two periodic clusters for the same values of the parameters 
r and e, but corresponding to different initial conditions. 

The inverse problem of determining the global coupling function in experimental systems is relevant since in general 
the specific functional form of the acting coupling is not known. This can be a complicated problem because, in 
addition, the exact form of the local dynamics may not be extracted in most situations. However, if the local 
dynamics is known some insight on the function H of a globally coupled system can be gained within the framework 
presented here. For example, in the case of a dynamical system showing two period-two clusters with partition p, the 
resulting asymptotic orbit [Hi,H2] can be obtained by measuring the cluster orbits and using Eqs. ( |l0|) and (|lj). 
For different realizations of partition p, the curve [Hi, H2] can be drawn as a function of p on the plane {Hi, H2), and 
compared with curves [Hi, H2] corresponding to known functional forms H such as those in Fig. 2. 

III. DRIVEN MAP ANALOGY 

As shown in Ref. ||ll[, the clustering behavior of GCM systems can be analyzed through its analogy with the 
dynamics of a single map subjected to an external drive. In order to interpret the results of the preceding section, let 
us consider an associated driven map 

St+i^{l-e)f{st) + eLt, (12) 

where st is the state of the map at discrete time t, f{st) is the same local dynamics as in Eq.(|l|), and Lt is an external 
driving term assumed to be periodic with period P. We denote the sequence of P values adopted by the periodic 
drive Lt by [Li,L2, ■ ■ ■ ,Lp]. The analogy between a GCM system and a driven map arises because in the former 
system (Eq. (p)) all the elements are affected by the global coupling function in exactly the same way at all times, 
and therefore the behavior of any element Xt{i) in the GCM is equivalent to the behavior of a single driven map (Eq. 
([l2|)) with Lt — H and initial condition Sq — Xoii)- Additionally, if a GCM system reaches a clustered, periodic 
collective state, its corresponding coupling function H follows in general a periodic motion. Thus the associated 
driven map (Eq. (|l2|)) with a periodic drive Lt should display a behavior similar to that of an element belonging 
to a periodic cluster in the GCM system. In particular, periodic drives resulting in periodic orbits of st in Eq. (03) 
may be employed to predict the emergence of clustered, periodic states in a GCM (Eq. (||)), regardless of the specific 
functional form of the global coupling H and without doing direct simulations on the entire GCM system. 

The driven map is multistable; i.e., there can exist several attractors for the same parameter values r and e. 
Specifically, for a given periodic drive [Li, L2, ■ ■ ■ , Lp], the map St may reach a number of distinct asymptotic responses 
Stij); {j — 1,2, ... , J), all with the same period, and depending on initial conditions. The orbits of st with period M 



can be expressed as a sequence of values [si(j),S2(j), • ■ ■ :Sm{J)]- The correspondence between a GCM system (Eq. 
(|l|)) in a state of K clusters with period P and its associated driven map (Eq. (p^)) can be established when M = P 
and J — K. 

Using this analogy, the main features in Fig. 2 can now be explained. In terms of a driven map subjected to a 
period-two drive Lt = [Li, L2] and the same parameters r = 1.7 and e = 0.2 as in the GCM systems considered in Fig. 
2, the bounded region on that figure contains the values of Li and L2 for which the driven map st just possesses two 
distinct asymptotic orbits (J = 2) of period two {M — 2), denoted by [si(l), S2(l)] and [si(2), S2(2)]. For values of Li 
and L2 outside this diamond shaped region, the map st may also reach a number of asymptotic periodic orbits but 
none with both J = 2 or ill = 2, at least on the interval —1 < Li,L2 < 1. Because of the analogy drawn above, all 
period-two motions [Hi, H2] of permutable coupling functions in GCM systems given by Eq. ([l]) with r = 1.7, e = 0.2 
and displaying two clusters in period two will fall on this bounded region of the plane {Hi,H2). Equivalently, a 
collective state of two clusters in period two can emerge in a GCM system only if its global coupling function has an 
orbit [Hi,H2] with values Hi = Li and H2 = L2 contained within the bounded region of Fig. 2. The boundaries 
of that region vary as the parameters r and e are changed. It is out of the scope of the present work to establish 
how the bounded region for two clusters in period two, as well as other regions for different clustered, periodic states, 
depend on parameters. However, it is worth noticing that this dependence can serve to characterize GCM systems 
with permutable coupling functions, and that this characterization can be obtained by the sole use of an associated 
driven map. 

The analogy between a GCM system with a given coupling function iJ in a two-cluster state and an associated 
driven map can be carried further by defining an associated coupling function [ pT| 

eH = i/(st(l),st(2),p); (13) 

that is, Qh is similar to the reduced two-cluster coupling function, such as Eqs.(o)-(P) but with the arguments Xi(l) 
and Xt(2) from the cluster trajectories replaced by the driven map orbits St(l) and St(2). The associated coupling 
function links the dynamics of clusters in a GCM system to the dynamics of a single associated driven map. The 
function Oh depends on the functional form of the coupling function H , on the partition {p, 1 — p} among the two 
clusters in the GCM, and on the orbits st(l) and st(2), which themselves are function of the period-two drive [Li, L2] 
and the parameters r and e. Thus for fixed r and e, and a given H, we have Oh = Qh{Li,L2,p). The main point 
is that, an equivalence between a GCM system Eq. ([l]) in a two-cluster, period-two state, and an associated driven 
map, Eq. dl3), with a period- two drive occurs when the following conditions are fulfilled 

9^(51(1), si(2),p)=ii, (14) 

eif(s2(l),S2(2),p)=i2. (15) 

Eqs. (|l4|)-(|l5|) constitute a set of two nonlinear equations for Li and L2, for a given p. The solution [L\, L2] of Eqs. 
([l4|)-(|l5|) predicts that the GCM possesses a state characterized by the coupling function motion [Hi = Ll,H2 = L2] 
and cluster orbits [xi(l),X2(l)] = [si(l), 52(1)], [xi(2), X2(2)] = [si(2), §2(2)]. The succession of solutions [Li,I/2]j 
as p varies, yields the curve corresponding to a given H in Fig. 2. Thus, each curve on the plane {Hi,H2) is 
parameterized by the fraction p. Moreover, there exist solutions [L|,L2] to Eqs. (p^)-(p^) only for an interval of p. 
Therefore, the curves in Fig. 2 can, in principle, be calculated a priori by using an associated driven map and just the 
specific functional form of H in each case. The range of possible cluster sizes, described by the values of the fraction 
p for which exist solutions to Eqs. (|lJ)-(Qq), can also be predicted by this method. Similarly, the regions of period 
two in Figs. 3(a)-3(d) can be obtained by varying the parameter e and calculating the interval of p for which Eqs. 
(||)-(|l5|) have solutions. 

IV. PREDICTION AND STABILITY OF CLUSTERS IN SYSTEMS WITH CONSTANT GLOBAL 

COUPLING 

Another simple clustered collective state in GCM systems occurs when the coupling function H remains constant 
in time, i.e., H = C. This behavior may take place in a GCM system with a permutable coupling function when 
K clusters, each having N / K elements and period K , are evolving with shifted phases in order to yield a constant 
value for H . That is, if the periodic orbits of K identical-size clusters are cyclically permuting in time, the resulting 
H becomes constant. For those collective states, the behavior of any of such clusters in the GCM system can be 
emulated by an associated driven map subjected to a constant forcing Lt = C 1l6|. In the case of a GCM displaying 
two equal size clusters in period two, this situation corresponds to the intersection of H with the diagonal in Fig. 2. 
The cluster orbits are then related as Xt(l) = [xi(l)i X2(l)] — [o-,b] and Xt(2) — [xi(2), X2(2)] — [b,a]. On the other 



hand, the associated driven map with Lt = C has a unique asymptotic period-two orbit St = [a, /3] on a range of C, 
where a and (3 are functions of C. The associated couphng function Qh also simphfies in such case. For the reduced, 
two-cluster couplings in Eqs (|6|)-(^, the corresponding associated coupling functions become 

(.) (a, /3) = ^ (« + /?) (16) 

0(/)(«,/3)-^[/(«) + /(/3)] (17) 

eA.(a,/3) = i|a-/3| (18) 

es(a,/3) = |ar/2|/?r/2. (19) 

Then Eqs. (p^)-(p^) with Li = L2 = C reduce to the single equation 

eH{a,l3) = C, (20) 

which can be seen as an equation for C, for given values of the parameters e and r. The solution C = C* of Eq. ( [2(]| ) 
provides a complete description of the GCM state since then a — a(C,), b = /3(C») and H = C*. Figure 4 shows the 
bifurcation diagram of Sf , Eq. (|l2[), as a function of the constant drive ii = C up to period two, with fixed parameters 
r — 2 and e = 0.4. The fixed point region in this diagram corresponds to one stationary cluster (i.e., a synchronized 
collective state) in the GCM, Eq. (|l|), with constant H. The period-two window corresponds to the values a and /? 
adopted by the driven map on this range of C. Once a{C) and P{C) are known from the bifurcation diagram, the 
function Qh{C) associated to any global coupling function in a GCM can be readily obtained, assuming that the 
GCM is in a state of two equal size clusters {p — 0.5), evolving out of phase with respect to each other. In Figure 4, 
the Qh functions associated to the four global couplings in Eqs.(||)-(||) with fixed e are shown as function of C As 
stated above, the solutions C* to Eq. (^0|) correspond to states in GCM systems with a coupling function reaching a 
stationary value H = C^,. Thus, the intersections of the Qh curves with the diagonal in Fig. 4 give all the possible 
states of GCM systems that maintain a constant H, either with one stationary cluster (if the intersection occurs 
on the fixed point window of the bifurcation diagram of St) or with two clusters in period two (if the intersection 
occurs on the period two window of the diagram) . Note that both the backward and the forward mean field couplings 
have the same two-cluster, period-two solution (x) = (/) = C* = 0.416, but these couplings have different functional 
dependence on the constant drive L^ = C. The coincidence of the couplings (x) and (/) for p = 0.5 was already 
seen in Fig. 2(b). Similarly, the geometric mean coupling H = x gives only one two-cluster, period-two solution at 
X — C^. = 0.305. In contrast, the dispersion global coupling, H = Ace, has three intersections with the diagonal: 
one corresponds to the synchronized stationary state in the associated GCM, with Aa; = C* = 0, and the other two 
correspond to different clustered states of the GCM, each consisting of two equal size clusters in period two, with 
Aa; = C, = 0.083 and Aa; — C^, — 0.25, respectively. All the states predicted by the intersections of the different 
Qh with the diagonal in Fig. 4, except one, are readily found in simulations on the corresponding GCM systems 
for appropriated initial conditions in each case. Actually, for a GCM with the coupling H = Ax, the predicted two- 
cluster, period-two state with Aa; = C* = 0.083 is unstable: it is never achieved in simulations on the GCM system, 
even when the initial conditions are chosen very close to that state. What is observed, instead, is the evolution of the 
GCM system towards either the stationary one cluster state with Aa; = or the two-cluster, period- two state with 
Aa; = 0.25. Thus, in addition to being predicted by the solutions of the equation Qh = C, the observed clustered 
states of a GCM displaying constant coupling must be stable, which implies some stability condition on the solutions. 
It can be shown (see Appendix) that for coupling functions satisfying J2i=i dH/dxi ~ 0, the condition dGn/dC > 1 
at the intersection with the diagonal implies that the corresponding solution is unstable. This is the case of the 
global coupling H = Ax. Note that the solution at C* = 0.083 is the only one for which ^U > 1 in Fig. 4 and 
therefore it is unstable, independently of the cluster fraction p. For the coupling H = (/), a stability analysis of states 
consisting of two or three clusters in period three has been performed by Shimada and Kikuchi |17]| . However, the 
simple criterium for instability ^1^ > 1 does not apply for H — (/). 

Constant coupling functions may also occur in GCM systems with different cluster sizes; that is the case of a GCM 
possessing dispersion global coupling H — Ax and displaying two clusters with any partition {p, 1 — p}, as seen in Fig. 
1(b). Figure 5 shows the associated function 8ax with fixed e as a function of the constant drive C for several values of 
the fraction p. There exist a critical fraction p^ — 0.25 bellow which only one solution corresponding to one stationary 
cluster, i.e., synchronization, can appear in the GCM system. Above this critical fraction, two states, each consisting 
of two clusters in period two, are additionally predicted by the solutions C = C, of Qax = C. These solutions emerge 
as a pair: one solution is always unstable since dQn / dC\c, > 1 there, and the other is the one two-cluster, period-two 
collective state with constant Ax that is actually observed in simulations on the corresponding GCM system. For the 
fraction pc = 0.25, there is a two-cluster, period-two solution Bax = C* = 0.125 that is marginally stable. 



V. CONCLUSIONS 

Most studies on GCM lattices and other globally coupled systems have assumed mean field coupling. However, 
other forms of global coupling may be relevant in some situations. We have analyzed, in a general framework, the 
clustering behavior in GCM systems subjected to permutable global coupling functions by considering the dynamics 
of the coupling functions. 

We have shown that different GCM systems can be represented by the orbits of their coupling functions on a common 
space. For simplicity, only collective states in GCM systems consisting of two clusters in period two were considered. 
We have shown that the functional form of the global coupling in a GCM system determines the periodicity of its 
motion and the possible distributions of elements among the clusters. The existence of a well defined interval of possible 
partitions among two clusters, out of which no clusters emerge in the system, has been observed experimentally [[l5| . 
In experimental or natural situations where clustering occurs, the specific functional form of the coupling is in general 
unknown. The present study may be useful to obtain insight into the acting global coupling function in practical 
situations. 

We have employed a previously introduced analogy between a GCM system and a single externally driven map |l|] 
in order to give a unified interpretation of the observed clustering behavior of the GCM systems considered in this 
article. A periodically driven map with local periodic windows can display multiple asymptotic periodic responses 
which are similar to cluster orbits in a GCM system with permutable H. This analogy implies that dynamical 
clustering can occur in any GCM system with a permutable coupling function and periodic windows in the local 
dynamics. The presence of windows of stable periodic orbits in the local map is essential for the emergence of clusters. 
In fact, no clustering is observed in a GCM system if the local maps do not have periodic windows fia]; what is 
observed instead is synchronization or nontrivial collective behavior, i.e., an ordered temporal evolution of statistical 
quantities coexisting with local chaos. 

The associated coupling function derived from the driven map analogy is particularly simple to use in the prediction 
of clustered states in GCM systems with two equal size clusters and exhibiting constant global coupling H = C. The 
associated coupling function Qh can be directly constructed from the bifurcation diagram of the steadily driven map. 
The cluster states are obtained from the solutions of Eq. (g^) and can be represented graphically in a simple way. 
Although Eq. ( pO[ ) has been used for the case of two clusters in period two, it can also be applied to find GCM states 
consisting of K equal size clusters in period K. In addition, the associated coupling function Qh carries information 
about the sta bility of the predicted two-cluster states. In particular, for the family of coupling function satisfying 



property (A9), the stability condition of clustered states in a GCM with constant iJ = C, is directly given by the 
slope (IQh /dC\c, ■ The example of a GCM system with dispersion coupling function H — Ax reveals that a constant 
coupling can also be maintained by clusters of different sizes. Our method based on Eq. ( po|) also predicts successfully 
the cluster states in these situations. 

The driven map analogy suggests that the emergence of clusters should be a common phenomenon which can 
be expected in various dynamical systems formed by globally interacting elements possessing stable periodic orbits 
on some parameter range of their individual dynamics. The examples presented here show that progress in the 
understanding of the collective behavior of globally coupled systems can be achieved by investigating their relation 
to a driven oscillator. 

ACKNOWLEDGMENT 

This work has been supported by Consejo de DesarroUo Cientifico, Humanistico y Tecnologico of the Universidad 
de Los Andes, Merida, Venezuela. 

APPENDIX: ON THE RELATION BETWEEN STABILITY AND THE ASSOCIATED COUPLED 

FUNCTION 

Consider a general GCM system with any global permutable coupling function. Suppose that the system reaches a 
state consisting of two clusters. Then the dynamics of the system reduces to two coupled maps Eqs. (p^-(pT|), i.e., 

Xt+i(l)= il-e)f{xt{l)) + eH{xt{l),Xtm,p)^ F(xt(l), Xt(2)), .,,^ 

Xt+i(2)= (l-6)/fe(2)) + eif(x*(l),X*(2),p)= G(xt(l),X*(2)). ^^^^ 

If the two clusters are in period- two orbits, the stability of this collective state in the GCM is given by the eigenvalues 
of the product of Jacobian matrices 
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If the two clusters move out of phase, the asymptotic state of the GCM can be described by the two orbits 

Xt{l) = [xi(l),X2(l)] = [a, 6] and xt(2) = [xi(2), X2(2)] = [5, a], which satisfy 



6= (l-e)/(a) + ei/(a,5,p), 
a= {l-e)f{b) + eH{a,b,p). 



For the local dynamics f{x) = 1 — rx"^, one gets 
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Xta)=a,Xt(2)=b 
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Xtil)=b:Xti2)=a 
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Xtil)=b.Xt(2)=a 



Consider now the dispersion coupling function, Eq. (^. This coupling belongs to the family of functions of N variables 

N 



H{xi,X2, 



,Xn) 



^/le(Xj - (x)), 



(A8) 



i=l 



where he is any even function of its argument. It can be straightforwardly shown that this family of functions possesses 
the property 



2—1 

Therefore, in a two cluster state, any H in this family of global coupling functions satisfies 

dH{xtil),Xti2)) , dH{xt{l),Xt{2)) 



dxtil) 



dxt{2) 



0. 



(A9) 



(AlO) 



If the two clusters evolve out of phase with respect to each othe r, a nd additionally the GCM has a coupling H with 
property (|A9D, then the two eigenvalues of the matrix J in Eq. (A2) become identical and their value is 



A = 2re{e - l){aHb + bHa) + 4r^ab{e - 1) 



(All) 



The stability criterion of this state is given by the modulus of the eigenvalue A; that is, |A| > 1 (|A| < 1) implies that 
the state is unstable (stable). The values a and b are, respectively, the values of a and /3 at the intersection of the 
function Qax with the diagonal in Fig. 4. 

Let us analyze the relationship between the eigenvalue A and the derivative dQAx/dC at the intersection points 
with the diagonal in Fig. 4 or Fig. 5. In general. 



deaeda dedf3_ 

dC '"d^dC^'dpdC' 



(A12) 



where 



-1 + iJ 

2r(e-l)' 



(A13) 



/3 



-l-R 

2r(e-l)' 



and 



Since 



then 



i? 



-3 - 4re2C + AreC + 4r - 8re + 4re2) 



2n1/2 



aa 9^9 1 ai? 


e 


dC dC 2r{e - 1) dC 


~R' 


dQ e (dQ dQ\ 




dC R\da dp) ' 





Since Qax has the same functional form as H — Ax, then Oax also satisfies property (AlO), that is 

dQAx dQAx 



da 



and therefore 



de 



Ax 



dC 



(3 ' 

2e dOAx 
R da ' 



(AM) 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 



Let C = C* be a value of C corresponding to the intersection of Qax with the diagonal in Fig. 4. Then a(C») = a 
and P{C^) — b; and 8ak(C*) — H{a, b), which gives 



dQAx 



dC 



R 



Using the fact that Ha — —Hi, from Eq. ( |A10| ), the eigenvalue A becomes 

A = 2re(e - l)Ha{b - a) + Ar^ab{e - 1). 



(A20) 



(A21) 



Eqs.(Al3 ) and (A14) with C — C* give the values a and b, respectively. Then, substitution of these values and Ha 



from Eq.(A20) in Eq. (A.21) yields 



A = i?^ 



dO/ 



dC 



(A22) 



Therefore, the condition dQ/dC\c, > 1 implies that |A| > 1, and thus the two-cluster, period-two solution with 
C* ~ 0.07 given by the intersection of Qax with the diagonal in Fig. 4 is unstable. Similarly, the solutions C* of 
©Ax = C for the different curves in Fig. 5 for which dQ/dC\c, > 1, are unstable. 

Note that the above stability result for H = Ax is also valid for any global coupling function satisfying property 



(A9) 
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FIG. 1. Dynamics of the coupling H = x (triangles) in a GCM system Eq. (|l|) with parameters r = 1.7, e — 0.2, displaying 
two clusters, each in period two. Cluster orbits are Xt(l) = [xi(l)iX2(l)] (circles) and Xt(2) = [xi(2), X2(2)] (squares), (a) For 
the partition p = pi = 0.4 and p2 ~ 0.6, H follows a period-two motion, adopting the values [Hi,H2]- (b) For {p = 0.5), the 
two clusters evolve out of phase with respect to each other and H remains constant at the value H = C = 0.3037. 
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FIG. 2. a) Curves of period-two orbits [Hi,H2] on the plane {Hi,H2) as p varies for the four coupling functions Eqs. 
(N)"® i'l corresponding GCM systems displaying two clusters. Parameters are the same for the four systems, r = 1.7, e = 0.2. 
The boundaries of the region where period-two orbits of any permutable H may take place is indicated with dashed lines, 
b) Magnification of a). Labels identify the curve associated to each H and the numbers besides the marks along each curve 
indicate the corresponding values of the cluster fraction p. The range of possible values of p for the different curves is also 
displayed on the figure. 
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FIG. 3. Regions of periodic motions on the plane (p, e) for different couplings H in GCM systems displaying two clusters. 
Local parameter is fixed at r = 1.7. The numbers on each region indicate the period of H on that region, a.) H = {x). b) 
H = (/). c) H = Ax; lA: there is only one stationary cluster with constant H — along the line p = 0, bistability occurs 
on the edge-shaped region: a state of one stationary cluster with H = coexist with a state of two out of phase clusters with 
constant, nonvanishing H; IB: there are two out of phase clusters giving constant H j^ 0. d) H — x. 
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FIG. 4. Bifurcation diagram of the driven map st, Eq. (|12|), with Lt = C, as a function of C. The asymptotic orbits of st 
are drawn with solid lines. The values a and /3 on the period-two window of the driven map are indicated; a = /3 on the fixed 
point window. The associated coupling functions Oh from Eqs. (|lq)-(|l9|) are also shown vs. C The intersections with the 
diagonal Qh ~ C are indicated by black dots and they correspond to the solutions C, of Eq. (EOl) 
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Intersections with the diagonal give the solutions C, . Solutions for which ^|c. > 1 are unstable. 
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